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Dissipative particle dynamics (DPD) is now a well-established method for simulating soft matter 
systems. However, its applicability was recently questioned because some investigations showed an 
upper coarse-graining limit that would prevent the applicability of the method to the whole meso- 
scopic range. This article aims to reestablish DPD as a truly scale-free method by analyzing the 
problems reported by other authors and by presenting a scaling scheme that allows one to gauge 
a DPD-simulation to any desired length scale. 
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I. INTRODUCTION 

Dissipative particle dynamics (DPD) was introduced 
in 1992 by Hoogerbrugge and Koelman^ as a novel 
method for performing mesoscopic simulations of com- 
plex fluids. Since then, the method has gained signif- 
icant theoretical support and refinement,^''*-^ and has 
been applied to fluid dynamics in numerous research 
areas such as rheology /'^''^'^ material sciences/^ilS and 
molecular biology, where membranes,^^'^^ vesiclesr^^ 
and micellar systems^^'^^'^^ have been modeled. 

Initially, DPD was understood to be a truly meso- 
scopic method able to bridge the whole gap between 
the underlying atomistic scale (in the range of nanome- 
ters and nanoseconds) that is accessible by molecular 
dynamics (MD) simulations and the macroscopic scale 
(in the range of micrometers and milliseconds) consid- 
ered by continuum descriptions. To fulfill this promise, 
it is crucial that the method is scalable, meaning that 
its coarse-graining level can be adjusted without intro- 
ducing serious artifacts that would render the method 
worthless. 

Coarse-grarnrng in DPD translates into having a num- 
ber V of physical molecules be represented by a single 
DPD particle.12 By N, we denote the total number of 
DPD particles in a simulation and it holds vN = iVphys, 
with A'phys the number of physical molecules the simu- 
lation represents. The main objective of this article is the 
comparison of DPD simulations with different coarse- 
graining levels V and v' . This motivates the introduction 
of the scaling ratio — N/N' ~ v' jv. In the follow- 
ing, functions of will be used to describe the scaling 
of various quantities at different coarse-graining levels. 
(One may set v = \ ■so that scaling ratio and coarse- 
graining level coincide and thereby simplify the nota- 
tion and providing a direct link to physical constants 
such as the Boltzmann constant; this is at the price of 
using DPD concepts at a length scale on which their in- 
terpretation becomes difficult.) 

It was originally stated that the method is scale-free, 
meaning that the parameters used in the simulation 



do not depend on the level of coarse-graining.'' In a 
later publication,^^ this earlier finding was declared er- 
roneous, and it was proposed that interaction parame- 
ters determining the conservative forces between DPD 
particles scale linearly with (\> when one goes from one 
coarse-graining level v to another v' . 

Based on this scaling relation, the performance 
of DPD was analyzed for various coarse-graining 
levelsji^^ It was found that there exists an upper coarse- 
graining level above which the simulated fluid freezes. 
Trovimof reported that this coarse-graining limit is dis- 
appointingly low,^^ and only allows up to about 10 wa- 
ter molecules to be grouped together into one DPD par- 
ticle. This limit would prevent DPD from covering the 
whole mesoscopic range and confines its applicability 
essentially to the order of magnitude of MD simula- 
tions. Dzwinel and Yuen even concluded that the DPD 
method would be best suited for the simulation of va- 
pors and gases (where the freezing artifact would hap- 
pen only for much higher coarsening levels). 

In this article, we argue that the original statement 
that DPD is scale free is in fact correct, and that the 
later refinement declaring that the DPD interaction pa- 
rameter scales linearly with the scaling ratio 4> is based 
on an inappropriate scaling assumption. More specifi- 
cally, we show that if a scaling method that preserves 
physical quantities is used and the DPD calculations are 
performed in reduced units (which is standard practice) 
then the results are independent of the level of coarse- 
graining. Expressed in more technical terms, if a phys- 
ically consistent scaling of all simulation parameters is 
maintained, the velocity increments calculated from in- 
tegrating the equations of motion of the DPD particles 
(c.f. Eqns. ^ below) do not depend on the coarse- 
graining value used for the DPD particles. This is in 
contrast to earlier publications. ^^'^^ 

In this paper, the physical molecules will refer to builk 
fluid particles (such as water), which is in accordance 
with the cited literature. Systems including surfaces, 
such as binary fluids, will be commented on in the dis- 
cussion. 
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FIG. 1: Starting from Frame A, a renormalization or scaling is 
performed that changes the interaction parameters, but does 
not change the units, leading to Frame B. This step also in- 
volves changing some of the force parameters as well a sys- 
tem constant (namely the Boltzmann factor) in order to main- 
tain the same behavior as the initial system had. The over- 
all system scaling is then changed with parameters being ex- 
pressed in terms of reduced units. The dynamics in the re- 
duced unit system (Frame C) is numerically exactly the same 
as in Frame A. In fact. Frame C is effectively a zoomed version 
of Frame A. 



The article is structured as follows: in Sec. Ill Ai we 
discuss the results of Groot and Rabone^^ and specify 
where we deviate from their analysis. We show that 
their approach of decreasing the number of DPD par- 
ticles (namely the particle density) while keeping rele- 
vant properties, in particular their radius of interaction, 
constant, is not appropriate. The alternative scaling pro- 
cess we employ is schematically shown in Fig. lA and 
IB. When we change the level of coarse-graining for the 
DPD-particles, we accordingly scale their number and 
adjusting their size (radius of interaction). Fig. IC de- 
picts the main result to be shown in this article, namely 
that by employing the correct scaling relations and unit 
reduction one has a complete equivalence of a simula- 
tion performed at the scale of Fig. lA with the zoomed 
version in Fig. IC. In other words, we demonstrate that 
one single DPD-simulation represents the dynamics of 
a whole family of physical systems at different levels of 
coarse-graining. 

This demonstration is split into two parts. In Sec. Ill Bl 
and III CI we change the interaction parameters accord- 
ing to the necessities of the scaling procedure we adopt. 
The interaction parameters have do be changed such 
that a system with many DPD-particles is mapped onto 
one with fewer, but larger and heavier particles. This 
is done without changing the physical properties of the 
complete system, namely the compressibility and ther- 
mal energy. Sec. Ill DI presents simulation results that cor- 
roborate the correctness of the derived scaling relations. 

In Sec. Uni we analyze the behavior of the DPD algo- 
rithm when the rescaled system is expressed in its own 
set of reduced units. The change of units affects the time 
scale and we will show that as a consequence, the nu- 
merical values that appear in the reduced unit system 



(Fig. [if) are identical to those of the original one and, as 
a consequence, their d5rnamics are equivalent. 

In order to present our argument, we need to define 
the nomenclature used and introduce some notation. By 
"coarse-graining" we understand the operation of coa- 
lescing I' physical particles into one DPD particle. By 
"scaling" we refer to the functional relation between 
the respective parameters of two systems with different 
coarse-graining resolutions. 

Conventionally, DPD operates in reduced units, such 
that energy is measured in units of kT {T being the tem- 
perature and k a scaled version of the Boltzmann con- 
stant to be introduced later), length in units of Tc, and 
mass in units of m, the mass of a single particle; in these 
units, length, mass, time and energy are dimensionless. 
We will use two different sets of reduced units, one for 
the coarse-graining resolution ly, and another one for i^'. 
We denote hy X a quantity expressed in reduced units 
with respect to coarse-graining level u. With X' we de- 
note the corresponding value of the same quantity cal- 
culated in the imits of X but with a different level of 
coarse-graining. In other words: moving from X to X' 
denotes the scaling operation. Finally, by X we refer to 
the quantity X' in its own set of reduced units given by 
the coarse-graining resolution ly'. That is, going from X' 
to X denotes the reduction of units. One may summa- 
rize the combined coarse-graining and change of units 
with the following diagram: 



X 



Scaling 
by factor t 



Reduction 
of units 



X 



Excellent descriptions of the DPD-method are given 
in various articles; we will not recapitulate the method 
itself, but instead refer to Groot and Warren.'' Here, we 
give only the definitions of the conservative, dissipative, 
and random forces in order to define the notation of the 
parameters: 



F,^ = -7c.^(n,)[(v.-v,) 



TT-R 



(1) 



where is the Euclidean distance between particles 
i and j, Vij is the unit vector pointing from particle j 
to particle i, and Xij equals one for pairs of particles 
separated by distances less than the force cut-off radius 
Tc and equals zero otherwise. The parameter de- 
termines the magnitude of the conservative interaction 
and will be regarded in this work as being the same for 
all pairs of particles: = o. is a random vari- 
able with Gaussian statistics, a vanishing mean and a 
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variance of 1/Ai for the numerical time step At (see 
e.g. Gardiner^^). As a consequence, the unit of Qj is 

time^^''^. Furthermore, the values of Qj in two different 
time intervals are uncorrelated. In general, if and 
{k, I) are different pairs of particles, dj and (ki are inde- 
pendent; however, in order to guarantee the centrality 
of all forces, one must require that Qj ~ (ji. 

The following dissipation-fluctuation relation^ leads 
to a thermal equilibrium at a given temperature T: 



Hr)? 



(2) 



Without loss of generality, we may take uj-^{t) — 
[u^{r)]^ for the dimensionless weighing functions 
uj^{r) and uj^'{r). The dissipation-fluctuation rela- 
tion (|2]l then reduces to 



(3) 



One is free to choose either uj-^{r) or uj-^'{r) without 
changing the thermodynamic equilibrium, but it is cus- 
tomary in the literature to take uj^{r) = 1 — r/vc mim- 
icking the conservative force F*^. One also notes that 
through Eqns. ^ and the relations between the pa- 
rameters therein are dependent on the temperature. 



II. SCALING DPD 

A. Compressibility and Equation of State 

Following Groot and Warren,^ we analyze the scaling 
behavior of the conservative interaction parameter a by 
relating the thermodynamic definition of the isothermal 
compressibility kt to the equation of state (involving a) 
of a system of DPD-particles. The isothermal compress- 
ibility is defined as the fractional change in the volume 
V that results from a change in the pressure P of the sys- 
tem, in a process where the temperature T is constant: 



}_dV_ 
VdP 



1 dp 

pap 



(4) 



Here p denotes the number density of particles, i.e. 
N/V. It is convenient to define the dimensionless pa- 
rameter by 



1 



1 dP 



ksTpKT ksT dp 
which we require to be invariant under scaling. 



(5) 



1 



ksTpKT 



1 dP 



ksT dp 



1 dP 



ksT dn 



const. 



exp 



(6) 

Here, n = vp denotes the molecular number density of 
the physical system and the subscript "exp" relates to 
the experimental value, while "sim" refers to the simu- 
lation value. 



The equation of state relates the pressure with the par- 
ticle number density p. From simulation results, Groot 
and Warren'' evaluated the virial expression 



P^pksT+^{y^n,F^X 



\J>1 

= pkeT+^p^ J^\a(^l~ g{rydr. (7) 

Here, g{r) denotes the radial distribution function. For 
densities p > 2 in reduced units (i.e. more than 2 par- 
ticles in a cubic box with linear dimensions of Vc), the 
following equation of state is a good approximation to 
the numerical simulations!^ 



P = pksT + aap 



(a = 0.101 ±0.001). (8) 



From this, one concludes that the part of the pressure 
caused by the conservative interaction scales linearly in 
a. From Eqn. we obtain that at constant temperature 



1 d 
k^dp 



{ksTp 



aap 



const. 



(9) 



Using p — n/v where n is the molecular number density 
and V the chosen coarse-graining parameter, it follows 
that 



2aan 
ksTv 



= const. 



(10) 



Since n and T are constant, Groot and Rabone^^ con- 
cluded that a must scale linearly with (p = v' /v in order 
to maintain isothermal compressibility under a change 
of the coarse-graining level. 

Several authors regard this scaling as an inherent 
drawback of the DPD-method, since on the microm- 
eter scale the method would appear to be effectively 
thwarted. We claim that this scaling - and hence its 
implications - is wrong. The problem is rooted in the 
construction of the equation of state (Eqn. |8]l. In their 
computer simulations, Groot and Rabone^^ decreased 
the density of DPD particles while keeping the inter- 
action cutoff radius Tc constant. This approach allows 
one to keep the system behavior invariant by scaling 
the interaction parameters shile changing the DPD par- 
ticle density, without altering their properties. How- 
ever, proper scaling means lowering the number of em- 
ployed DPD particles while simultaneously enlarging 
their interaction radius. The difference is illustrated in 
Fig. 121 FramelJJ'V is taken to be a system with fine coarse- 
graining. Frame [2^ represents a scaled system with a 
lower DPD particle density but unchanged particle di- 
ameters. The result is that the mutual overlap of the soft 
particles is smaller (as seen in Frame|2^). Hence it is in- 
tuitively clear that the interaction parameter has to be 
increased in order to keep the system properties con- 
stant; formally, this argument is reflected in Eqn. i fTOl l. 
In contrast. Frame llfZ shows the system with the same 
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scaling ratio as for Frame |2^, but with the relative over- 
lap of the interacting particles kept constant, which is 
accomplished by scaling j'c- A closer examination of 
Frame |2{Z! shows that it is part of a magnified version 
of Frame |2j\, namely a system where all the lengths as- 
sociated with a single DPD particle have been uniformly 
scaled by a factor (p while keeping the system size con- 
stant {L' ~ L). This results in the following scaling re- 
lations for the coarse-graining level, number, mass, and 
size of DPD particles: 

v' = (jyv 
N' = (t)-^N 
m! = (j)m 

= ^c, (11) 
where d is the number of dimensions of the system (see 

Fig.Ei::). 



B. Scaling of the Potential Energy 

We start by calculating the change of potential energy 
t/ of a system of DPD particles enclosed in a box that 
undergoes compression. This change is related to the 
compressibility of the system and is required to be in- 
variant under scaling. In practice, we require the de- 
pendence of a on (f), such that the chosen coarse-graining 
level does not affect the compressibility. For the uncom- 
pressed system, we have 

^o = Ef^(----^)'- (12) 

This equation holds for soft core repulsions which are 
used throughout the literature. In general, the poten- 
tial may be viewed as a harmonic approximation of any 
potential close to an energy minimum. For an isotropi- 
cally compressed system with box length (1 — e)L, where 
e <C 1 is the relative compression parameter, the change 
in the interparticle distance Ar^ (e) is not assumed to be 
the same for all pairs of particles. However, we require 
that 




FIG. 2: Schematic of the scaling process: Frame A shows a 
DPD simulation with a cutoff radius of Tc- Frame B depicts the 
coarse-graining procedure performed in Groot and Rabone^^. 
With changing particle density, the particle diameter is kept 
constant while the interparticle force is increased to maintain 
the system pressure. However, in order to properly conserve 
systemic parameters like compressibility, both the interaction 
parameter and the interaction cut-off radius need to be in- 
creased as measured in physical units. Frame C depicts the 
proper scaling in coarse-graining. Along with a decrease of 
the particle density (in physical units), the interaction range is 
increased. In this case, the interaction parameter a scales dif- 
ferently than in Frame B in order to preserve systemic proper- 
ties. 



Ary(e) = ery +0(e2), 



(13) 



which means that we rule out (first order) phase transi- 
tions under compression. The total potential energy of 
the compressed system is then given by 



(14) 



To first order in e, we obtain for the change of internal 
energy 



Af/ = C/,-f/o = Ex^,a(l 



(15) 



Because the change rn internal energy of the system as a 
whole has to be invariant under scaling, we have 



^Xy O I 1 - ^ 



N' 

E 



1-^1 er'. 



(16) 



Due to the scaling of A^, the number of terms in the sum 
of the left hand side of Eqn. ((T6] | is proportional to 0^^. 
Since we require AC/ to be invariant imder scaling, the 
force constant a has to scale as 



(17) 



the <\)^ coming from the change in the number of terms 
in the sum and the (p^^l'^ from the change in length 
scale. We note that this scaling differs from that found 
by Groot and Rabone where a scaled linearly with (\) (see 
Eqn. l|6ll). This scaling is the result of our requirement 
of maintaining the fractional particle overlap during the 
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change of the coarse-graining level. Everything else be- 
ing equal, this scaling would seemingly still imply an 
upper coarse-grarnrng limit, although not as severe as 
the initial result of Groot and Rabone^^. However, scal- 
ing affects not only length scales and the interaction pa- 
rameter, but also the energy and time scales implicit in 
the simulation. 

C. Scaling of the Kinetic Energy 

Since the unprimed and the primed systems should 
be physically equivalent, we further require that T' = T 
and Q' = Q, where Q is the thermal energy, i.e. heat 
content, of the whole system; with these requirements, 
physical properties become independent of the coarse- 
graining, and the coarse-graining level is solely a simu- 
lation parameter. In order to be consistent with statis- 
tical mechanics, a proper consideration of the reduction 
of the number of degrees of freedom resulting from co- 
alescing several physical particles into a DPD particle is 
required. The heat content in a 3-D system consisting of 
particles with no internal structure is given by 

Qphys = = ^(™p,,,V^,^J, (18) 

where /phys — dNp^ys is the number of degrees of free- 
dom. Forming N DPD-p articles, i.e. coherently moving 
groups of I' physical particles, changes the relation be- 
tween kinetic energy, temperature and number of parti- 
cles. It holds 

Q-{fcr=^(mv2) (19) 

with / = 3iV, k = vks, N = -/Vphys/i^, rn — z^mp^ys. Con- 
sequently, from the required invariance Q' = Q = Qphys 
and T' — T for two simulations, we get a scaling relation 

k' = (t>k (20) 

We introduce the parameter k that has the role of a Boltz- 
mann constant for systems with reduced number of de- 
grees of freedom. (This scaling affects the dissipation- 
fluctuation relation Eqn. ^ which now reads ct^ = 
2kTj) That Eqn. | [T9l l and hence Q = Q' is consistent 
with the DPD method is discussed and explicitly cor- 
roborated in Sec. Ill DI 

The behavior of 7 and a imder scaling remains to be 
determined. Examining Eqn. |[7|, we note that the pres- 
sure is independent of 7 and tr. This means that, with re- 
spect to static compressibility, we have significant free- 
dom in the choice of the scaling function F: 

7' = r(0)7 

a' = (r(0)</.)i/V (21) 

the latter equation is a consequence of a' = ^/l-^'k'T' . 
Dimensional analysis motivates the choice 



which in turn implies that 

7' ^ 0i-i/''7 

a' = ^i-^/'^a. (23) 

This specific choice will later be shown to be crucial for 
establishing the scalability of the method, but investiga- 
tions with other goals (such that deciding on a coarse- 
graining level V with least artefacts) may require alter- 
native gauges. 

To summarize what has been established through- 
out the last two sections, changing the level of coarse- 
graining in DPD requires the following scaling relations: 

m' = (j)m 1' — 4>^'^^'^^ fr..s 

= 4>ydr, a' = ^i-^/'^a ^ ' 

k' = (j)k 

D. Experimental validation 

In order to illustrate the correctness of the above scal- 
ing arguments, we have measured the pressure in sim- 
ulations as a function of a for different coarse-graining 
levels and temperatures. In these simulations, the DPD 
particles are confined to a box with hard walls. When a 
particle collides with a wall, it is reflected elastically and 
the instantaneous impulse normal to the wall is mea- 
sured. The pressure is measured as the time-average of 
the normal forces on the walls divided by the surface 
area of the cube: P — {niAv±/ (AAt)) where A is the 
area of the box and Av_l is the component of the parti- 
cle velocity orthogonal to the wall. 

In Fig.|3]we show a comparison between rc = 1 (white 
squares), corresponding to normal DPD, and a rescaled 
system with r'^ = 2rc, corresponding to = 8 diamonds 
with black dots). Particle density and box length are 
set to Sr^'^, lOrc, and ir'^'^, 10?'^, respectively. Different 
temperatures Ti, T2, and T3 are considered, such that 
kTi = 1, kT2 = 2, fcTg = 3, and k'T{ = 1, k'T^ = 2, 
fc'Tg = 3, respectively (bottom to top in Fig.|3); this is a 
slight deviation from the conventional usage of reduced 
units. For the rescaled simulation, all parameters {N, m, 
rc, a, k, 7, and a) have been scaled according to Eqn. l l24l l. 
We want to show that the pressure of this system is in- 
variant under scaling, for all values of a and T, if we 
follow the described scaling relations. 

In Fig. |3] the pressure is plotted as a function of 
a/0^~^/'', so that we should obtain the same curve inde- 
pendent of the coarse-graining level. On the one hand, 
for a = 0, the case of an ideal gas, the pressure should 
be given by 

P\a=o = r^Nc^kT/V = pkT, (25) 

which it indeed is seen to be. On the other hand, kinetic 
gas theory establishes 



r(0) = </.i-i/'i, 



(22) 



P = p(tov2)/3. 



(26) 
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This constitutes a relationship between the thermostat 
and the conservative mechanical interactions with the 
walls, which is non-trivial for our simulation. The 
way we measure the pressure gives a direct relation 
to (v^). The pressure calculated from kinetic quanti- 
ties (Eqn. | |26|| ) coincides with the value required from 
the thermodynamic relation Eqn. ( [191 . This shows that 
Eqn. (19)1 holds for our simulation and thereby confirms 
the scaling relation for k given in Eqn. l|20l l. 

Eqn. (|8) estimates a linear relation between a and P. 
This relation is shown as solid lines, whereby the Boltz- 
mann factor in Eqn. l|8ll has been replaced by its scaled 
version k. For each line, ^l^^g is the theoretical value 
and the slope is the average obtained from each data set 
(P(a) — ^'(O)) /a. Unless a is too large, the pressure in- 
creases linearly with a. For high values of a, one finds 
that Eqn. l|8ll is not a good approximation to the mea- 
sured pressure, anymore. 



III. REDUCTION OF UNITS 

Our goal is now to show that the velocity increments 
Av obtained from integrating the forces are unchanged 
when the scaling is combined with the according reduc- 
tion of units: Av = Av', which implies that the relative 
particle motions are unaffected by scaling in the reduced 
unit systems. 

In the DPD literature, length, mass, and energy are 
considered primary units, leading to a derived unit of 
time r given by 




From the above arguments we obtain the time-scale in 
the rescaled system: 



The plot also shows that the scaling of a is correct as 
the pressure values of the original and the rescaled sys- 
tems fall right on top of each other, regardless of the 
value of a. The conclusion from these comparisons is 
that the pressure is indeed invariant when we use the 
scaling relations given in Eqn. | |24|| . 



k'T 



' (j)m 



(/.^/''t. (28) 



50 



Thus, length and time scale in the same way. This also 
implies that the velocities are invariant as they are given 
by Ar/At. 

The random variable Qj has the unit t^^/^, as noted 
in the discussion following Eqn. [T] Given the scaling of 
T, it follows that 




FIG. 3: The calculated pressure as a function of the interaction 
strength a. We show the pressure for the basic DPD model, 
corresponding to rc = 1 (white squares), and a rescaled sim- 
ulation of the same 3D system with 4> = 8 corresponding to 
r'c = 2rc (diamonds with black dots). Particle density and box 
length are set to Sr^T^, lOrc, and Sr^"^, lOr'^, respectively. Dif- 
ferent temperatures Ti, T2, and T3 are considered, such that 
kTi = 1, kT2 = 2, fcTs = 3, and k'Tl = 1, k'T^ = 2, fc'Tg = 3, 
respectively (bottom to top). Solid lines show the prediction of 
the equation of state (Eqn. ^) with a estimated from the sim- 
ulation results (for the case of r-c — 1). The pressure increases 
approximately linearly for a not too large. Furthermore, the 
values for the original and the rescaled simulations coincide. 



and therefore 



I/-' 

ij 



-1/(2^)/^.. 



(29) 



(30) 



Since the other terms in Eqn. ^ are all scale-free, 
the three force components of DPD all scale by a factor 

When velocity increments are calculated during one 
time step, one finds that the force scaling is canceled by 
the scaling of mass and time. This is shown below for 
the conservative force: 



< ^ m' mn' < ^ \ 



m 

(j) m 
= AvF. 



(31) 



The calculations for [Avf^]' and [Av^*]' give similar re- 
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stilts: 



(j) m 



= Avf (32) 

= 2^ ^ (^uOC^ji-^j At 

(33) 



= Avf 



Since Af = Ar'/r^ and At = At'/r', we get for the ve- 
locity increment by considering 



Av, = Af j/At = Av' r'/^ 



(34) 



Because time and length scale in the same way we get 
r'/r^ = t/tc- Combining this with Av = Av' (Eqns. (|3Tl l 
to | [33l l) one finally obtains 



which implies 



Av, = Avi 



m = r(t). 



(35) 



(36) 



What remains to be shown is the scaling of the re- 
duced parameters a, 7, and a. Since a scales like en- 
ergy over length, when going to the reduced units of the 
primed system, we have 



Ll-1/dj.l/d 



k'T' 



a-r-= = a, 



(j) kT 



(37) 



and similarly, since 7 scales like energy over length and 
velocity, from 7 = 7' r'^ / {k'T'r') we get 



Ll-l/dj.2/d J2 



^ ^ k'T'r' 



7^7^ = 7. (38) 



From the fluctuation-dissipation relation it follows 
again that 



a = a. 



(39) 



Hence, scaling and unit reduction precisely cancel each 
other. As a result, the DPD formalism is scale-free. This 
means that the calculation with one and same set of pa- 
rameter values represents systems at arbitrary lengths 
scales. 

Note that in order to understand the according phys- 
ical time scales we have to comment on transport prop- 
erties of the method and the scaling of fluctuations. 
We base our argument on diffusion, but could equally 
well consider viscosity since the two are related by the 



Schmidt-number, which is dimensionless. The numer- 
ical equivalence of the measured diffusion constants 
D = D (being a consequence of the equality of the mea- 
sured displacement) and the fact that diffusion scales 
like length squared over time causes an apparent prob- 
lem: it seemingly implies that relative fluctuations D/L 
stay constant instead of vanishing. This problem dis- 
appears when one calibrates the simulation to an actual 
physical system. Assume that the cutoff radius is related 
to a physical length by ~ I [cm] . We then have 



cm' 



phys 



sec 



(40) 



with D 



phys 



^^^^'^^ ^ referring to the diffusion con- 



stant in physical units. We get 

D 



phys 



/'[sec]. 



(41) 



Consequently, expressed in physical units, it holds for 
the fluctuations 



lim 



D, 



phys 



^phys 



1 



(42) 



IV. SUMMARY AND DISCUSSION 

We have shown that each individual term of a re- 
duced unit DPD calculation is scale-free, and hence so is 
DPD as a whole. We emphasize that this result was only 
achieved by a combined effect of scaling, the reduction 
of units, and finally by requiring a specific but obvious 
scaling of 7 and a. 

The scaling of energy per particle (or equivalently the 
number of degrees of freedom in the system) is shown 
to require increasing the value of k by the same factor 
of (f> that the number of particles is reduced by. This 
scaling of k is required to keep the total system energy 
constant. The increase in k compensates for the freezing 
effect seen by other authors who failed to take this into 
account and who also did not uniformly scale all length 
parameters. 

In this article, we have shown that the DPD-method 
is scale-free for the simulation of bulk fluids. This is not 
necessarily the case e.g. for binary mixtures of liquids 
A and B where several conservative interaction param- 
eters occur, e.g. aAA,ciAB,ciBB- Whereas bxilk interac- 
tions given by aAA, clbb scale as discussed in this arti- 
cle, aAB is a surface term that determines interf acial en- 
ergy and therefore scales differently. This is, however, 
beyond the scope of this article. 
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